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Climate impacts of ice nucleation 

A. Gettelman, 1 ’ 2 X. Liu, 3 D. Barahona, 4 ’ 5 U. Lohmann, 2 and C. Chen 1 

Received 16 April 2012; revised 11 September 2012; accepted 16 September 2012; published 19 October 2012. 

[i] Several different ice nucleation parameterizations in two different General Circulation 
Models (GCMs) are used to understand the effects of ice nucleation on the mean climate 
state, and the Aerosol Indirect Effects (AIE) of cirrus clouds on climate. Simulations 
have a range of ice microphysical states that are consistent with the spread of observations, 
but many simulations have higher present-day ice crystal number concentrations than 
in-situ observations. These different states result from different parameterizations of ice 
cloud nucleation processes, and feature different balances of homogeneous and 
heterogeneous nucleation. Black carbon aerosols have a small (—0.06 Wm~ 2 ) and not 
statistically significant AIE when included as ice nuclei, for nucleation efficiencies within 
the range of laboratory measurements. Indirect effects of anthropogenic aerosols on 
cirrus clouds occur as a consequence of increasing anthropogenic sulfur emissions with 
different mechanisms important in different models. In one model this is due to increases in 
homogeneous nucleation fraction, and in the other due to increases in heterogeneous 
nucleation with coated dust. The magnitude of the effect is the same however. 

The resulting ice AIE does not seem strongly dependent on the balance between 
homogeneous and heterogeneous ice nucleation. Regional effects can reach several 
Wm . Indirect effects are slightly larger for those states with less homogeneous 
nucleation and lower ice number concentration in the base state. The total ice AIE 
is estimated at 0.27 ±0.10 Wm~ 2 (1 cr uncertainty). This represents a 20% offset of the 
simulated total shortwave AIE for ice and liquid clouds of —1.6 Wm~ 2 . 

Citation: Gettelman, A., X. Liu, D. Barahona, U. Lohmann, and C. Chen (2012), Climate impacts of ice nucleation, J. Geophys. 
Res., 117, D20201, doi:10.1029/2012JD017950. 


1. Introduction 

[ 2 ] Ice clouds play an important role in the climate system. 
Ice clouds reflect solar radiation back to space, cooling the 
planet. However, ice clouds, being cold, radiate much less 
long wave radiation to space than low clouds or the surface, 
thus also heating the planet. The balance of these two 
processes determines the net ice cloud radiative forcing, 
which is finely balanced between warming and cooling, with 
warming thought to be slightly larger. 

[ 3 ] Changes to ice cloud microphysics might alter this 
balance. Since ice clouds are nucleated by aerosol particles, 
changes to the aerosol composition in ice cloud regions may 
alter cloud microphysics and ice cloud radiative forcing. 
These are commonly called ‘Aerosol Indirect Effects’, or 
AIE [Twomey, 1977], reviewed by Lohmann and Feichter 
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[2005]. While the AIE of liquid clouds is predominately a 
short wave cooling by reflection of solar radiation off of 
brighter cloud tops, cold ice clouds also have longwave 
effects, and complex interactions between aerosol and crys- 
tal number \Kdrcher et al., 2006] that can lead to warming or 
cooling. There remain major uncertainties in ice nucleation 
processes [DeMott et al., 2011], and while several repre- 
sentations of ice nucleation now exist [Liu and Penner, 
2005; Kdrcher et al., 2006; Barahona and Nenes, 2009], 
and global models are incorporating them with comprehen- 
sive aerosol schemes [Lohmann and Hoose, 2009; Gettelman 
et al., 2010] there are limited estimates of the global impact 
of ice nucleation on climate. Penner et al. [2009] performed 
estimates using off-line radiative transfer calculations and 
found no net effect from sulfate, but a cooling effect from 
Black Carbon soot (see below). Liu et al. [2009] estimated 
the ice AIE using the Community Atmosphere Model (CAM) 
version 3 of 0.5-0. 7 Wm -2 and lower stratospheric water 
vapor increases of 0.3 ppmv. 

[ 4 ] There are major uncertainties in the types of ice nuclei 
and the mechanisms by which ice clouds form that com- 
plicate assessment of ice indirect effects. The role of dif- 
ferent aerosol types in nucleating ice at different threshold 
humidities is uncertain [DeMott et al., 2011], Attempts at 
estimating indirect effects need to account for these uncer- 
tainties, which are driven by wide variations in ice micro- 
physical measurements [Kramer et al., 2009]. 
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Table 1. Description of Simulation Pairs Used in This Study a 


Name 

Ice Nucleation 

Description 

CAM5-LP 

Liu et al. [2007] 

Base CAM5 (LP2007 Ice Nucleation) 

CAM5-FixedIN 

Cooper [1986] 

Ice Nuclei a function of temperature 

CAM5-LP-BC-0. 1% 

Liu et al. [2007] 

Base + 0.1% BC included as heterogeneous IN 

CAM5-LP-BC-2% 

Liu et al. [2007] 

Base + 2% BC included as heterogeneous IN 

CAM5-LP-BC- 1 00% 

Liu et al. [2007] 

Base + all BC included as heterogeneous IN 

CAM5-BN 

Barahona and Nenes [2009] 

BN2009 with D 504 > 0.05 pm 

CAM5-BN-Low 

Barahona and Nenes [2009] 

BN2009 with D S04 >0.1 pm 

CAM5-LP-Low 

Liu et al. [2007] 

Base with 40% less ice mass 

CAM5-LP-W10 

Liu et al. [2007] 

CAM5-LP2007 with sub-grid W x 10 

CAM5-LP-FixedDN 

Liu et al. [2007] 

Base code with fixed drop number 

CAM5-BN-FixedDN 

Barahona and Nenes [2009] 

BN2009 code with fixed drop number 

ECHAM5-FixedDN 

Kdrcher and Lohmann [2002] 

Fixed drop number 


a Each simulation pair is the difference between aerosol emissions from the year 2000 and 1850, with annually repeating climatological SST and 
greenhouse gases. For CAM5-FixedIN, Cooper [1986] specifies IN = 209 L _1 for T < — 30°C. 


[5] One of the major uncertainties with respect to ice 
nucleation is understanding the importance of natural and 
anthropogenic Black Carbon (BC, also known as soot). 
Process [Kdrcher and Lohmann, 2003; Kdrcher et al., 2007] 
and global [. Lohmann et al., 2004; Hendricks et al., 2005; 
Penner et al., 2009; Liu et al., 2009; Hendricks et al., 201 1] 
model studies found significant impacts of BC on cirrus 
clouds if they are assumed to be efficient ice nuclei [e.g., 
Penner et al., 2009]. Barahona et al. [2010] compared sev- 
eral descriptions of heterogeneous ice nucleation within the 
same parameterization framework and found that a soot ice 
activation fraction of 1% would be enough to significantly 
impact the global distribution of ice crystal concentration 
at cirrus levels. 

[6] The assumption that BC can heterogeneously nucleate 
ice crystals at humidities below the homogeneous nucleation 
limit (RH, of 130-150% depending on temperature) is not 
certain. Cloud chamber measurements by Mohler et al. 
[2005] found that 0. 1-0.3% of uncoated BC particles 
freeze at supersaturations from 1.1- 1.3 (110% < RH,- < 
130%). But DeMott et al. [2009] could not find heteroge- 
neous freezing of biomass burning aerosols in controlled 
burning laboratory experiments. Neither did Friedman et al. 
[2011] for controlled experiments with coated soot at pure 
ice formation temperatures (— 35°C). However, assuming 
that BC aerosols are efficient ice nuclei at low super- 
saturations (RH, = 130%), Penner et al. [2009] find using 
off-line calculations significant and large impacts of 
anthropogenic BC on cirrus clouds, with radiative forcing of 
up to —0.4 Wm“ 2 from the indirect effects of BC, but with 
a high BC mass loading of 3x observations [Wang et al., 
2009]. Crawford et al. [2011] found that between 0.1-1% 
of soot particles with low organic carbon content would 
nucleate ice at T ~ 220 K in the deposition mode. Barahona 
[2011] showed that this is consistent with a soot ice 
nucleation fraction increasing from about 0.01% to 10% for 
RH, ~ 1 35% to RH,- ~ 1 70%, respectively, and with the soot 
IN activity increasing for decreasing temperature. 

[7] This study will investigate (a) the different climate 
states of humidity and clouds simulated by different ice 
nucleation schemes, (b) anthropogenic effects on ice clouds 
and (c) the role of BC in altering ice cloud radiative forcing. 

[8] We will investigate these properties using two different 
advanced General Circulation Models (GCMs), the National 
Center for Atmospheric Research (NCAR) Community 


Atmosphere Model version 5 (CAM5), and the European 
Center Hamburg (ECHAM) model version 5 with the 
Hamburg Aerosol Model (ECHAM5-HAM). We will use 
3 different ice nucleation parameterizations in several dif- 
ferent configurations, and a version of CAM5 with no ice 
indirect effects, to examine the impact of ice nucleating 
aerosols on climate. 

[ 9 ] Models and experiment methodology (Section 2) 
are defined below. Baseline results with present-day emis- 
sions, and the impact of BC in present day are discussed 
in Section 3. Ice AIE estimates are presented in Section 4. 
The paper finishes with a discussion (Section 5) and con- 
clusions in Section 6. 

2. Models and Experiments 

2.1. CAM5 Description 

[ 10 ] The primary model used in this study is CAM5 from 
the National Center for Atmospheric Research (NCAR). 
CAM5 includes an advanced physical parameterization suite 
[Gettelman et al, 2010; Neale et al., 2010] that is well suited 
for understanding ice nucleation. CAM5 has a 2-moment 
cloud microphysics scheme [Morrison etal., 2005; Gettelman 
et al. , 2008], coupled to a Modal Aerosol Model with 3 modes 
[Liu et al, 2012a]. The model allows ice supersaturation, and 
links ice cloud particles consistently to aerosols through ice 
nucleation (see below). Aerosols in CAM5 do not interact 
with convective cloud drops and ice crystals. Convective 
detrainment of ice into anvils is based on temperature (all ice 
below — 30°C) and a prescribed mean radius of 25 / im . 
CAM5 has a consistent treatment of the radiative effects of 
ice clouds and snow (see Gettelman et al. [2010] for details). 

2.1.1. CAM Ice Nucleation 

[ 11 ] CAM5 simulates ice nucleation in mixed phase and 
cirrus (pure ice) clouds. Ice nucleation for cirrus (pure ice) 
clouds (T< — 35°C) is parameterized following Liu and 
Penner [2005] and Liu et al. [2007], hereafter LP2007, 
implemented as described by Gettelman et al. [2010]. This 
simulation we term CAM5-LP (Table 1). The parameteriza- 
tion features competition between heterogeneous nucleation 
(immersion freezing) on dust with homogeneous nucleation 
on sulfate. Unlike LP2007, nucleation of BC is not active in 
the standard CAM5 [Gettelman et al., 2010]. CAM5 homo- 
geneous ice nucleation uses Aiken mode sulfate aerosol 
particles larger than 0.1 /im diameter. CAM5 uses dust in 
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the coarse mode as potential heterogeneous ice nuclei (IN) 
[Gettelman et al., 2010], consistent with classical theory and 
IN observations [Pruppacher and Klett, 1997]. We also test 
adding BC as IN with the same properties of dust as in Liu 
et al. [2007] or Penner et al. [2009], with an efficiency of 
0.1%, 2% or 100% (CAM5-LP-BC cases in Table 1). CAM5 
in the mixed phase regime (— 40°C < T < — 3°C), features 
empirical formulations for ice nucleation for deposition and 
condensation freezing [Meyers et al. , 1992], immersion 
freezing [Bigg, 1953], prescribed without explicit dependance 
on aerosols, and contact freezing on dust [Young, 1974]. 

[ 12 ] Since the grid-scale vertical velocity is the net of 
larger small scale motions, CAM5 ice nucleation is driven 
by a sub-grid scale vertical velocity. The velocity is derived 
from the Turbulent Kinetic Energy (TKE) scheme [Park and 
Bretherton, 2009]. In the upper troposphere, this velocity is 
typically 5-10 cm s -1 , smaller near the tropopause, and 
larger near deep convective clouds. In general this velocity is 
lower than that seen in limited observations from aircraft in 
cirrus anvils [Jensen et al., 2009] of 30-50 cms' 1 , and an 
order of magnitude smaller than convective motions, but is 
reflective of mean motions in large scale cirrus clouds and 
radiative uplift [Cord et al., 2005], We have performed a 
sensitivity experiment with the sub-grid vertical velocity for 
ice nucleation increased by a factor of 10 (CAM5-LP-W10), 
where the mean is 30-80 cm s -1 to test potential impacts 
of the specification of vertical motion on the results. 

2.1.2. Alternative Ice Nucleation 

[ 13 ] For comparison, we have also implemented an alter- 
native cirrus ice nucleation scheme to the base version of 
CAM5. This is the parameterization of Barahona and Nenes 
[2009], hereafter BN2009 (CAM5-BN). Similar to LP2007, 
it estimates the maximum supersaturation achieved in a rising 
air parcel and an aerosol spectrum to estimate ice crystal 
concentrations. However, whereas LP2007 is derived from 
fits to the numerical solution of the parcel model equations, 
BN2009 is derived from an approximation to their analytical 
solution. Using offline calculations Barahona and Nenes 
[2008] have shown that for similar values of the ice deposi- 
tion coefficient and pure homogeneous freezing, ice number 
concentration (TV,-) from BN2009 and LP2007 agree within 
a factor of 2 for T > 200 K. However, the two parameteriza- 
tions largely differ in their treatment of heterogeneous ice 
nucleation and the competition between homogeneous and 
heterogeneous freezing. LP2007 used a combination of clas- 
sical nucleation theory and empirical formulations to describe 
the freezing of dust and soot IN. BN2009 uses a generalized 
ice nucleation spectrum, which can have any functional form. 
Moreover, the maximum ice nuclei concentration that allows 
homogeneous freezing (used to scale TV,- as a function of the 
IN concentration in both parameterizations) in LP2007 is a 
function of the vertical velocity only, whereas for BN2009 it 
also depends on the characteristics of the ice nucleation 
spectrum, the dust and BC concentrations, and the cloud 
formation conditions. Thus BN2009 in principle allows a 
greater flexibility in describing heterogeneous ice nucleation. 
In this work the heterogeneous ice nucleation spectrum 
from Phillips et al. [2008] is used to describe ice nucleation 
on soot and dust IN in the deposition, condensation and 
immersion modes. We have adjusted the dust surface area in 
the Phillips et al. [2008] spectrum following Eidhammer 
et al. [2010]. The scheme use Aiken mode sulfate aerosol 


particles larger than 0.05 /jm diameter and dust in the coarse 
mode as potential ice nuclei as with LP2007. The sulfate size 
threshold has been altered from LP2007 in CAM5-BN to 
ensure reasonable ice number concentrations and consistent 
global cloud radiative forcing. 

2.2. ECHAM5-HAM Description 

[ 14 ] For comparison to CAM5, we also perform simulations 
of ice nucleation using a similar methodology with a different 
GCM. ECHAM5-HAM is a GCM with similar capabilities 
to CAM5. It has a 2-moment cloud microphysics scheme 
with prognostic number and mass for liquid and ice 
[Lohmann et al., 2007] coupled to a 7 mode aerosol model 
(HAM) [Stier et al., 2005]. Aerosol effects on ice clouds are 
parameterized as detailed in Lohmann et al. [2008] using the 
ice nucleation parameterization of Kdrcher and Lohmann 
[2002], If the dust number concentration exceeds 10 L -1 
heterogeneous nucleation sets in above a critical relative 
humidity of 130%. Otherwise cirrus form by homogeneous 
nucleation. BC is treated as an immersion IN in mixed phase 
clouds in ECHAM-HAM. 

2.3. A1E Analysis 

[ 15 ] There are several different ways to analyze aerosol 
indirect effects. To isolate effects for ice is difficult, as 
aerosols can also be active as cloud condensation nuclei 
(CCN), and also may have direct effects (scattering radiation 
back to space, or absorbing it in the case of BC). So we will 
use several different techniques to estimate the indirect 
effects in the simulations. Some of the complexities are 
described by Lohmann et al. [2009]. 

2.3.1. RFP Methods 

[ 16 ] The most straightforward method diagnoses a total 
aerosol effect by the Radiative Flux Perturbation (RFP) in 
the Top Of Atmosphere (TOA) radiative flux (R) between 
a simulation with present-day (year 2000) emissions of 
aerosol and aerosol precursors and one with pre-industrial 
(1850) emissions (RFP = R 2000 — Risso)- Both runs have the 
same Sea Surface Temperature (SST) and Greenhouse Gas 
(GHG) climatology (representing the years 1980-2000). The 
aerosol indirect effect can also be equated with a change in 
cloud radiative forcing (CF = R a// — R (/ ,. for LW and SW) 
between two runs (AIE = dCF LW + dCY S w) where R r/j . is the 
TOA flux calculated without clouds and R a// is the flux for 
all conditions. In this work we will use ‘c/’ to indicate a 
change between 2000 and 1 850 with the same code (within 
one experiment), and ‘A’ to indicate differences between 
different experiment pairs. To isolate the effect of ice nucle- 
ation, we can difference simulation pairs with ice nucleation 
linked to aerosols (e.g., CAM5-LP) and simulations in which 
ice nucleation is estimated only as a function of temperature 
(CAM5-FixedIN). The FixedIN simulation follows Cooper 
[1986] as in Gettelman et al. [2008], This can be estimated 
with the RFP so that AIE,- ce = ARFP = RFP^p — RFF^edW- 
Since the direct scattering effects of aerosols may also change 
between simulations, it is necessary to factor out changes 
to the clear sky shortwave flux (dFSC) that may vary between 
pairs of simulations (AFSC), so that A I E,- ce = ARFP — 
AFSC. AIE can also be estimated with the changes in 
cloud radiative forcing (CF) so AIE,- ce = ACF = c/GF LP — 
dCF FixedIN . These 3 estimates are shown in Table 5. 
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2.3.2. Fixed Drop Number 

[ 17 ] Another approach is to fix the liquid drop number 
(eliminating liquid indirect effects) and just focus on the ice 
phase. We have also done this in the simulations, performing 
several experiments with fixed drop number (Table 4, 
FixedDN). The AIE of ice is then t/RFP = dR = R 2 ooo — R 1 xso, 
minus the clear sky effect of aerosols, estimated by the change 
in clear sky shortwave flux (c/FSC), or AIE /ce = c/RFP — c/FSC. 
The AIE for a fixed drop number run can also be estimated by 
just the change in CF between 2000 and 1850 (AlE, ce = <7CF). 

[is] These different methods provide a spread of estimates 
in AIE. Fortunately, we will show that they are broadly 
comparable and provide a consistent set of numbers. 

2.4. Experiment Description 

[ 19 ] Each experiment listed in Table 1 and discussed 
below is two simulations, with different values of aerosol 
emissions representative of the years 2000 and 1 850. All other 
inputs and boundary conditions are kept constant at year 2000 
conditions (such as SST and GHGs). CAM5 simulations have 
horizontal resolution of 1.9° x 2.5° and are run for 6 years, 
analyzing the last 5 years. This is sufficient to get robust 
global statistics, but regional changes may not be signifi- 
cant. CAM5 simulations are performed using LP2007 ice 
nucleation (CAM5-LP), IN as fixed function of temperature 
(C AM5 -F ixedIN), BC added to LP2007 with different effi- 
ciencies (CAM5-LP-BC cases), and BN2009 ice nucleation 
with two different size thresholds for sulfate particles 
(CAM5-BN and CAM5-BN-Low). We also perform simu- 
lations with an alternative ice tuning for LP2007 (CAM5- 
LP-Low) that has similar global cloud forcing but 40% less 
ice to accommodate changes to liquid clouds. A sensitivity 
test was performed with sub-grid vertical velocity for ice 
nucleation increased by a factor of 10 (CAM5-LP-W10). 
We also perform fixed drop number experiments with 
LP2007 (CAM5-LP-FixedDN) and BN2009 (CAM5-BN- 
FixedDN), and with ECHAM5-HAM (ECHAM5 -F ixedDN ) . 
Drop number is not allowed to vary from a prescribed value 
set to produce reasonably present-day cloud forcing. 

3. Results: Ice Nucleation 

3.1. Evaluation of CAM 

[ 20 ] CAM5 has been extensively evaluated for the ice phase 
in the upper troposphere against observations by Gettelman 
et al. [2010]. The CAM5.1 model used here performs essen- 
tially identically with respect to ice mass and ice supersatu- 
ration. The simulations for the base case are very similar to 
Liu et al. [2012b], with the exception that Liu et al. [2012b] 
increase the sulfate (SO4) in ice nucleation by including all 
sulfate, while here and in Gettelman et al. [2010] we limit 
sulfate in ice nucleation to D > 0. 1 pm in the CAM5-LP case. 
We will show a few examples of ice cloud microphysics from 
the simulations, and reference other work that has validated 
the base case climatology. 

[ 21 ] Here we look briefly at the climatology of cirrus cloud 
microphysics in the base CAM5 (CAM5-LP) simulations. 
Figure 1 illustrates annual mean ice effective radius (Figures la 
and lb) and number concentration (Figures lc and Id) for 
the upper troposphere at 2 levels. 142 hPa (Figures 1 a and 1 c) 
features tropical cirrus clouds in the Tropical Tropopause 
Layer (TTL), while 232 hPa (Figures lb and Id) also includes 


cirrus clouds in middle and high latitudes. For tropical 
cirrus at 142 hPa (Figures la and lc), simulated Nj averages 
1 10 L _1 (Table 2) with a range from 100-300 L -1 . Simulated 
effective radii (Figure la) are from 20-50 pm. At lower 
altitudes and warmer temperatures (232 hPa), simulated 
ice Nj is generally lower, 100-200 L _I over land and only 
50-70 L _I over ocean (Figure Id). Simulated size is larger, 
and peaks at 70 pm in the tropics, and 45 pm in midlatitudes. 
Ni in these simulations is lower than in Liu et al. [2012b], 
due to the difference in sulfate number used. 

[ 22 ] To compare these simulated values to observations, 
Figure 2 illustrates the number concentration (Figure 2a) 
and (number weighted) effective radius (Figure 2b) as a 
function of temperature. Collected aircraft observations for 
tropical and midlatitude cirrus reported by Kramer et al. 
[2009, Figure 9] are also shown in Figure 2 (gray symbols) 
with error bars indicating 25th through 75th percentile, and 
‘X’ indicating the median (50th percentile). Percentiles fol- 
low Wang and Penner [2010], and effective radii are area 
weighted to compare to the model (they are multiplied from 
what appears in Kramer et al. [2009] by a factor of 3 to 
account for difference from number weighting assuming an 
exponential distribution). Also shown are observations of 
thin tropopause level cirrus from Lawson et al. [2008] (black 
square) and tropical convective outflow cirrus from Jensen 
et al. [2009] (black triangle). Simulated values are monthly 
means. In most cases, the most frequently observed value 
is the 25th percentile (lower concentration in Figure 2a), 
as noted in Salzmann et al. [2010], CAM5-LP (orange) Nj 
(Figure 2a) is slightly higher than observed at low tempera- 
tures, and lower at higher temperatures than observed. 
This indicates a different gradient than observations, likely 
due to the predominance of homogenous nucleation at 
lower temperatures. Some of the simulations (e.g., CAM5- 
BN-Low) do have Nj closer to observations (50-150 L” 1 at 
1 90 K). As a result, the ice effective radius (number weighted), 
is generally smaller than observed at low temperatures, but 
similar to observed at higher temperatures. Above 210 K, 
agreement between CAM5-LP and observations is good, 
with a tendency toward smaller radius and higher Nj at 
low temperatures. This might be due to suppression of 
homogenous freezing [Jensen et al, 2005]. ECHAM5-HAM 
(red) has higher Nj and smaller sizes at warmer temperatures. 
The spread of simulations allows a test of the impact of 
these differences on the resulting ice indirect effects. 

[ 23 ] As illustrated by Gettelman et al. [2010, Figure 10], 
CAM5 generates reasonable values and frequency of occur- 
rence of ice supersaturation when compared with available 
in-situ aircraft or satellite observations. Maximum values of 
about 150% RHi are seen, which is close to the threshold 
for homogeneous nucleation (150-160% at the coldest tem- 
peratures). These values are similar to Liu et al. [2012b], 
but the different balance of homogeneous and heteroge- 
neous nucleation alters some of the details at high ice 
supersaturation. 

[ 24 ] As shown by Liu et al. [2012b] and Gettelman et al. 
[2010], CAM5 generally is able to reproduce important radia- 
tive properties of clouds. Cloud radiative forcing and cloud 
fraction compare favorably to satellite observations [Gettelman 
et al., 2010, Figure 12; Liu et al., 2012b, Figure 8], as do ice 
mass in the atmosphere, and the partitioning between liquid 
and ice phases with temperature [Gettelman et al., 2010, 


4 of 14 



D20201 


GETTELMAN ET AL.: ICE NUCLEATION 


D20201 


A)CAM5-LP 142 hPa R e Ice B)CAM5-LP 232 hPa R e Ice 



C)CAM5-LP 142 hPa Num Ice D)CAM5-LP 232 hPa Num Ice 



Figure 1. Maps of annual mean (a and b) simulated ice effective radius (R e ) and (c and d) ice number 
concentration at 232 (Figures lb and Id) & 142 (Figures la and lc) hPa from the Base CAM5-LP 
simulation (year 2000). 


Figure 3]. These results indicate that the simulations are 
capable of representing the base state of cirrus and mixed 
phase clouds. 

[ 25 ] We performed a test to understand the sensitivity of 
the results to the assumed sub-grid vertical velocity used to 
drive ice nucleation (CAM5-LP-W10). Results indicate a 
near doubling of ice crystal numbers (Table 2) to values that 
are significantly higher than observations, and several W m~ 2 
increase in the gross cloud forcing (brighter clouds). 

3.2. Impact of Ice Nucleation 

[ 26 ] Now we look at the changes to the simulations from 
different ice nucleation approaches. We first look at the 


change between the CAM5-FixedIN case and CAM5-LP 
case that incorporates the effects of aerosols. We then 
explore differences with the CAM5-BN case. We look also 
at the impact of adding BC to the ice nucleation in CAM5- 
LP-BC, as originally included in LP2007, but removed in 
CAM5 [Gettelman et al., 2010] and the CAM5-LP case. 

[ 27 ] Figure 3 compares the grid box averaged ice (not 
including snow) mass (Figures 3a and 3b) and in-cloud ice 
number (Figures 3c and 3d) between CAM5-LP (Figures 3a 
and 3 c) and the ice nucleation fixed as a function of temper- 
ature following Cooper [1986], CAM5-FixedlN (Figures 3b 
and 3d). All else being equal in the simulations, the CAM5- 
FixedlN simulation has less ice mass and significantly lower 


Table 2. Key Cloud Properties in Base Case (Year 2000 Runs) a 


Run 

SWCF (Gbl) 
(WnW 2 ) 

LWCF (Gbl) 
(Wnr) 

IWP (Gbl) 
(g m ~ 2 ) 

INC 150 
(L ‘) 

REI150 

(t im ) 

Observed 

-47.1 

26.5 




CAM5-LP 

-51.8 

24.0 

17.7 

111 

31 

CAM5-FixedIN 

-49.2 

18.5 

13.1 

11 

41 

CAM5-LP-BC 0.1% 

-51.5 

23.6 

17.5 

102 

31 

CAM5-LP-BC 2% 

-57.8 

31.5 

22.3 

184 

27 

CAM5-LP-BC 100% 

-78.3 

51.8 

33.9 

2180 

20 

CAM5-BN 

-51.6 

23.7 

16.9 

252 

27 

CAM5-LP-Low 

-51.7 

23.5 

10.7 

131 

29 

CAM5-BN-Low 

-49.0 

20.7 

15.5 

34 

40 

CAM5-LP-W10 

-54.7 

27.7 

20.4 

200 

29 

CAM5-LP-FixedDN 

-53.7 

24.1 

10.8 

135 

28 

CAM5-BN-FixedDN 

-53.1 

27.1 

11.4 

686 

24 

ECHAM5-FixedDN 

-52.5 

27.3 

4.8 

986 

32 


a Shown are global Short Wave Cloud Forcing (SWCF), Longwave Cloud Forcing (LWCF) and Ice Water Path (IWP). Also shown are tropical 
(30S-30N) averages at 150 hPa of ice effective radius (REI150) and ice number concentration (INC 150). Observations of cloud forcing are from updated 
version 2. 1 of the CERES-EBAF data described by Loch et al. [2009], 
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Figure 2. Annual mean simulated ice (a) number concentration (TV,) and (b) area weighted effective radius 
(. R e ) from year 2000 simulations of different experiments (colored lines). Error bars indicate 1 standard 
deviation of individual points in a temperature bin. Also shown are observations from Kramer et al. 
[2009] (gray) with error bars indicating the 25th and 75th percentile range, with median (50th percentile) 
indicated by ‘X’. Also shown are estimates in sub-visible cirrus by Lawson et al. [2008] (squares), and 
observations in tropical anvil clouds by Jensen et al. [2009] (triangles). 


ice number in the tropical upper troposphere at cirrus altitudes 
(300-100 hPa). The tropical average TV, at 142 hPa in CAM5- 
FixedlN is one order of magnitude lower than the CAM5-LP 
case (Table 2). These differences are not unexpected since 
the CAM5-Fixed IN run has a constant activated IN concen- 
tration of 209 L -1 for T < — 30°C, generally less than either 
CAM5-LP or CAM5-BN. There are also significant differ- 
ences at high latitudes, with less ice in polar regions. In the 
Northern Hemisphere, the reduced ice in the CAM5-FixedIN 
simulation has higher liquid mass below it (not shown). In the 
CAM5-FixedIN case, reduced ice sedimentation and precip- 
itation as snow collects less liquid at lower levels and allows 
liquid concentrations to be higher than the base case. This has 
significant impacts on the surface radiation budget in the 
Arctic (not shown). 

3.3. Effect of Alternative Ice Nucleation 

[ 28 ] There are also significant impacts when the BN2009 
scheme is used to nucleate ice heterogeneously and homo- 
geneously (CAM5-BN). A key feature of these simulations 
is that the BN2009 scheme has more active ice nucleation 
with higher ice numbers (2.5 x) than the CAM5-LP case 
(Table 2). This is because the threshold diameter for sulfate 
aerosols has been decreased from 0.1 (CAM5-BN-Low) to 
0.05 /itm (CAM5-BN) to maintain the energy balance at 
the top of the atmosphere. With the larger size threshold 
(the CAM5-BN-Low case), the BN2009 scheme produces 
lower ice numbers (more similar to observations by Kramer 
et al. [2009, Figure 2] and has very low LW cloud forcing 
when compared to observations from the Clouds and the 


Earths Radiant Energy System (CERES) satellite [Loeb et al. , 
2009], as illustrated in Table 2 (CAM5-BN-Low). We will 
examine the AIE of this low ice number case however, since 
at cold temperatures CAM5-BN-Low has ice numbers and 
sizes closer to observations (Figure 2). 

[ 29 ] With the BN2009 scheme (CAM5-BN), the ice 
nucleation has changed substantially. Figure 4 shows the 
fraction of activated ice nuclei from homogeneous freezing 
in the CAM5-LP case (Figure 4a) and the CAM5-BN case 
(Figure 4b). In CAM5-LP there is a balance between homo- 
geneous and heterogeneous nucleation. Particularly in the 
Northern Hemisphere, with high dust loading, there is sig- 
nificant heterogeneous nucleation in the upper troposphere. 
With the BN2009 parameterization as used in CAM5-BN, 
almost all ice nucleation in the upper troposphere is due to 
homogeneous freezing (Figure 4b). A similar result has been 
noted by Liu et al. [2012b]. The result is a climate simulation 
with similar radiative transfer and cloud forcing (Table 2), but 
a very different microphysical balance of clouds. 

3.4. Impact of Black Carbon 

[ 30 ] We also explore adding Black Carbon (BC) as an ice 
nucleus in the LP2007 scheme. This is done with three cases 
(CAM5-BC) at different efficiencies for BC as an efficient 
heterogeneous ice nucleus. Table 2 illustrates that there is a 
strong sensitivity to the fraction of BC that can serve as IN. 
Here we explore the assumption that either all (100%) or 
small fractions (2%, 0.1%) of BC are efficient ice nuclei, 
behaving like dust and nucleating heterogeneously at 
RH, = 130% as in LP2007. Even with 0.1% BC, there is a 
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Figure 3. Zonal mean (a and b) CAM5 simulated ice mass and (c and d) number concentration from the 
Base (Figures 3a and 3c) and Fixed IN (Figures 3b and 3d) year 2000 simulations. 


significant shift in the balance between homogeneous and 
heterogeneous nucleation from the base CAM5-LP case 
(Figure 4a) toward heterogeneous nucleation dominating, 
especially in the Northern Hemisphere (Figure 5a). With 2% 
BC available for activation, almost all ice nucleation occurs 
heterogeneously (Figure 5b). Table 2 indicates that for the 
case of 100% BC and 2% BC the gross cloud forcing is 
unrealistically high compared to observations from CERES 
observations [Loeb et al., 2009] shown in Table 2 by a factor 
of 1.5-2 (LW and SW) for the 100% case, and 1. 3-1.1 (LW 
and SW) for the 2% BC case, and a doubling (2%) or tripling 
(100%) of the ice mass, with large increases in ice number. 
For the 0.1% BC case (CAM5-LP-BC 0.1%), there is little 
impact on the simulated climate, and the energy balance of 


this run remains similar to observations and the CAM5-LP 
case. We also explored including BC ice nucleation in the 
BN2009 scheme (not shown). BN2009 assumes a reduced 
capability of BC as an IN, and also has higher homogeneous 
nucleation, so there was little impact seen in the simulation 
for adding BC as IN. 

4. Results: Aerosol Indirect Effects 

[ 31 ] Given this set of simulations, we now have several 
different model configurations with very different micro- 
physical balances, but similar climate states. Based on the 
spread in current observations (Figure 2) these states are not 
distinguishable from each other. Even the large differences 



Figure 4. Zonal Mean fraction of activated IN due to homogeneous nucleation in the year 2000 
(a) CAM5-LP and (b) CAM5-BN simulations. 
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Figure 5. Zonal Mean fraction of activated IN due to homogeneous nucleation in the year 2000 CAM5- 
LP-BC simulations with (a) 0.1% BC and (b) 2% BC efficiencies. 


between the BN and Base case ice number in the upper 
troposphere do not have a large impact on cloud radiative 
forcing (Table 2). Cloud radiative forcing is similar because 
the microphysical balance of clouds (ice number, habit and 
mass) combine in a non-unique function to generate similar 
radiative forcing from different states. The question remains: 
will these different states and different balances of ice 
nucleation result in differences in how ice clouds change, for 
example, in response to aerosol perturbations? 

[ 32 ] We use the CAM5 and ECHAM5-HAM experiments to 
test this hypothesis. We perform an additional set of experi- 
ments where the emissions of aerosols reflect emissions in 
1850. All other aspects of these simulations (GHGs and SSTs) 
remain the same. First we analyze the quantitative AIE of 
all clouds in the simulations, then discuss differences between 
the various simulations for ice AIE. The basic statistics of the 
differences between the 2000 and 1850 emission runs are 
shown in Table 3. The total liquid and ice AIE from the 
CAM5-LP case is 1.0 Wm" 2 , when estimated either from 
the change in cloud radiative forcing (cfCF = —1.11 Wm -2 ), 
or the change in RFP, adjusted for the direct effect of aerosol 
scattering (dR - c/FSC = -1.36 + 0.37 = -0.99 Wm -2 ). 

4.1. Difference From Fixed IN 

[ 33 ] First we examine the effect of adding ice nucleation 
from the FixedlN simulation (CAM5-FixedlN) where aerosols 
do not affect cirrus clouds. Adding ice nucleation (CAM5-LP) 


reduces the change in top of atmosphere flux (dR) and 
change in cloud forcing (CF) from the fixed IN case by 
+0.2 Wm -2 (Table 3). The pattern of the change is illustrated 
in Figure 6. The impact is mostly confined to the northern 
hemisphere (where aerosol emissions change), and there are 
significant differences over the N. Pacific and over the Arctic. 
High variability in the Arctic reduces significance of regional 
changes there in these short (5 year) simulations, but high 
latitude Arctic changes do not contribute much to the global 
value. Most of the global effect is due not to a change in ice 
mass, but to a moderation of the change in liquid mass 
(Table 3). Based on budgets and process rates, Ice nucleation 
appears to have a ‘glaciation’ effect, where formation of ice 
at higher altitudes results in snow fall and accretion of liquid 
by the falling snow. This is particularly important for the 
Arctic and mixed phase clouds. 

[ 34 ] In the base case simulation, the change of radiative 
forcing is due to a change in the ice nucleation regime. 
Figure 7 shows a map of the fraction of ice crystals nucleated 
by homogeneous freezing at 232 hPa from the CAM5-LP 
case for the year 2000 (Figure 7a) and 1850 (Figure 7b). 
Figure 7a corresponds to the same field as in Figure 4a. 
In 1850 homogeneous nucleation fractions are low in regions 
dominates by dust, especially over Europe, the Sahara, 
and Asia, extending into the N. Pacific. In 2000, many of 
these regions have significantly more homogeneous nucle- 
ation due to higher sulfate aerosol loading. In the TTL (140- 


Table 3. Radiative Property Changes (Year 2000-1850) From Simulations 3 


Run 

dR 

(Wm -2 ) 

dCF 

(Wm -2 ) 

r/LWCF 

(Wm -2 ) 

c/SWCF 

(Wm -2 ) 

c/FSC 

(Wm -2 ) 

c/IWP 

(gm -2 ) 

c/LWP 

(gm -2 ) 

c/INC 

(L -1 ) 

c/REI 

(/mi) 

CAM5-FixedIN 

-1.58 

-1.33 

0.64 

-1.97 

-0.29 

0.07 

4.7 

0 

0.1 

CAM5-LP 

-1.36 

-1.11 

0.54 

-1.65 

-0.37 

0.05 

3.1 

19 

-0.8 

CAM5-LP-BC 0.1% 

-1.37 

-1.25 

0.47 

-1.72 

-0.34 

0.13 

3.6 

14 

-1.3 

CAM5-LP-BC 2% 

-1.42 

-1.44 

2.45 

-3.89 

-0.29 

0.88 

4.8 

58 

-1.3 

CAM5-LP-BC 100% 

-1.96 

-0.80 

-1.71 

0.91 

1.60 

-1.0 

-0.4 

-100 

0.4 

CAM5-BN 

-1.31 

-1.25 

0.80 

-2.05 

-0.29 

0.26 

3.6 

60 

-1.5 

CAM5-LP-Low 

-1.25 

-1.16 

0.63 

-1.79 

-0.40 

0.18 

3.8 

26 

-1.2 

CAM5-BN-Low 

-1.25 

-1.10 

0.31 

-1.41 

-0.42 

0.02 

3.1 

5 

-1.6 

CAM5-LP-W10 

-1.02 

-0.81 

1.77 

-2.58 

-0.28 

1.14 

3.5 

45 

-1.3 

CAM5-LP-FixedDN 

-0.06 

0.31 

0.41 

-0.10 

-0.34 

0.05 

0.2 

20 

-0.9 

CAM5-BN-FixedDN 

0.10 

0.19 

1.16 

-0.97 

-0.20 

0.32 

0.8 

131 

0.0 

ECHAM5-FixedDN 

0.13 

0.39 

0.22 

0.17 

-0.32 

0.03 

0.14 

133 

-0.7 


“Illustrated are change in top of atmosphere radiative fluxes (R), net cloud forcing (CF) as well as the long-wave (LWCF) and shortwave (SWCF) 
components, the change in clear-sky shortwave radiation (FSC), ice water path (IWP) and liquid water path (LWP). Also shown are tropical (30S-30N) 
averages at 1 50 hPa of changes in ice effective radius (REI) and ice number concentration (INC). 
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Figure 6. AIE from ice clouds estimated by change in Top of Atmosphere (TOA) radiation (dR) between 
2000 and 1850 from (a) CAM5-LP and (b) CAM5-FixedIN simulations. 


100 hPa) ice number increases by 20% between 1850 and 
2000 (Table 3). Thus the balance of ice nucleation is shifting 
due mostly to increased upper tropospheric sulfate enhancing 
homogeneous ice nucleation. This affects ice number and 
radiative properties of cirrus clouds. 

[ 35 ] For small nucleation efficiencies (0.1%), the effect of 
BC on AIE is not large (Table 3, case CAM5-LP-BC 0.1%). 
For higher efficiencies (CAM5-LP-BC 2% or 100%) the base 
state of the simulations indicates significantly higher gross 
cloud radiative forcing (Table 2). Analysis of these runs for 
AIE is hence difficult. However, the runs feature significantly 
more heterogeneous nucleation due to soot freezing in the 
immersion mode. For high soot activation, this causes a 
qualitative change in the effect of anthropogenic emissions. 
CAM5-LP-BC 100% has reductions (not increases) in crystal 
number, and increases in crystal size (Table 3) due to the 
anthropogenic effects of aerosols, changing the sign of the 
perturbation to the LW and SW cloud forcing. This is 
probably because homogeneous freezing is limited due to 
more efficient heterogeneous freezing on BC. However, 
neither the 2% or 100% case has a ‘realistic’ climate. The 
definition of ‘realistic’ has many dimensions, but here we 
illustrate it by comparison to the overall cloud radiative 
forcing (CF) from observations in Table 2. CF is simply a 
fundamental indicator of the cloud thickness and micro- 
physics and strong gross forcing indicates clouds which 
differ significantly from observations (also seen in number 
concentrations in Table 2). 


4.2. Fixed Drop Number 

[36] Another way to isolate the impact of ice nucleation on 
clouds is to perform simulations where the liquid clouds have 
a fixed drop number, and droplet activation does not depend 
on aerosols. This eliminates any direct liquid cloud response 
to changing aerosol emissions (though semi-direct effects 
of direct aerosol heating may occur). Three simulations 
with fixed drop number were conducted to estimate AIE 
from ice clouds: CAM5 with LP2007 (CAM5-LP-FixedDN), 
CAM with BN2009 (C AM5 -BN-F ixedDN), and one with 
ECHAM5-HAM (ECH AM5 -F ixedDN). Table 3 indicates 
almost no change in top of atmosphere radiation in the 
FixedDN runs, because the SW cooling effects of aerosols in 
the clear sky (dFSC in Table 3) are also acting in these 
simulations to offset the change in cloud forcing. Hence we 
look at the net cloud forcing as a measure of the ice only AIE 
for these runs in Figure 8, using both LP and BN para- 
meterizations. In both cases in Figure 8 the patterns are noisy 
due to a small signal. However, most of the signal is in the 
N. Hemisphere with regional perturbations of 3 Wm~ 2 or 
more. For the CAM5-LP-FixedDN simulations, the signal is 
concentrated over Europe, E. North America, E. Asia and 
the N. Pacific, all regions of high anthropogenic sulfate 
emissions. The signal in the cloud forcing is a net warming 
of -0.3 Wm" 2 . 

[ 37 ] The change in cloud forcing between 1 850 and 2000 in 
a fixed drop number run has also been analyzed in ECHAM5- 
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Figure 7. Map of 232 hPa fraction of ice crystals nucleated due to homogeneous freezing in the CAM5 
Base simulations for (a) 2000 and (b) 1850. 
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B) CAM5-BN-FixedDN ACF 



Figure 8. AIE from ice clouds estimated from the change in Cloud Forcing (ACF) from CAM5 fixed 
drop number runs between 2000 and 1850 for (a) CAM5-LP-FixedDN and (b) CAM5-BN-FixedDN. 


FixedDN (Table 3). The change in cloud forcing is slightly 
higher than in CAM5-LP-FixedDN (0.4 Wm" 2 , Table 3), but 
the mechanism is different. In ECFIAM5-FixedDN, there is 
an increase in longwave forcing (dLWCF, Table 3) but also a 
decrease in shortwave forcing (positive dSWCF in Table 3), 
associated with higher ice crystal numbers and slightly 
smaller sizes (Table 3). This arises despite an increase in 
heterogeneous nucleation in present-day, because more dust 
coated by sulfate and soluble organics is available. In the 
LP2007 scheme (CAM5-LP, with or without BC), heteroge- 
neous immersion freezing occurs for coarse mode dust, 
independent of sulfate coating. Thus a change in the hetero- 
geneous immersion freezing of sulfate coated dust is not 
included in CAM5, so heterogeneous freezing of dust acts in 
the same way in CAM5 in 2000 and 1850. The ice water path 
is lower in ECHAM5-FixedDN (4.8 g m -2 ) than in CAM5 
(17.7 g m -2 , Table 2). An increase in ice crystal number in 
ECHAM5-FixedDN has a larger effect in the LW than the SW 
for thin clouds [Fusina et al., 2007]. However, in ECHAM5- 
FixedDN, the cloud cover is also slightly reduced (not shown) 
so that the shortwave cloud forcing decreases slightly (posi- 
tive dSWCF in Table 3). The longwave cloud forcing 
increases as the increase in ice water path with more numer- 
ous and smaller ice crystals outweighs the decrease in cloud 
cover. There is virtually no change in cloud cover in CAM5, 
so the effects are only due to cloud microphysical changes. 
There are thus two differences between CAM5 and ECHAM- 
HAM. First the treatment of BC as heterogeneous immersion 
freezing nuclei in ECHAM-HAM. Second, the net effect on 
cloud forcing in each model is the same, but for different 
reasons: thinner clouds are changing their LW and SW bal- 
ance differently in ECHAM5-HAM than in thicker CAM 
clouds, combined with effects of changes in cloud area cov- 
erage in ECHAM5-HAM. 

4.3. Net AIE 

[38] There are several different ways to calculate the net 
AIE of ice clouds in the simulations, as detailed in Section 2.3. 
The two basic methods are to estimate it from the fixed 
drop number simulations, where liquid clouds do not change, 
or to estimate the impact as a difference from the effects 
in the Fixed-IN run using ice nucleation as a function of 
temperature. The results are generally consistent between 
methods. We also explore whether significant differences 
exist in the AIE of cirrus clouds with the different ice 


nucleation parameterizations and microphysical balance, or 
with the inclusion of Black Carbon as IN. 

[ 39 ] First we estimate the ice indirect effect from the fixed 
drop number runs. This can be done in two ways, using 
either the change in cloud forcing (dCF) or the change in 
flux (dR) adjusted for the direct scattering effects of aerosols 
(dR - dFSC). Table 4 illustrates these values. For the LP2007 
simulation with fixed drop number (CAM5-LP-FixedDN), 
both values are essentially the same (~0.3 Wm -2 ). There is 
slightly more spread in the two values from the CAM5-BN- 
FixedDN simulation, but they both indicate a wanning 
due to anthropogenic aerosols of 0.2-0. 3 Wm -2 . The gross 
change in cloud forcing in the BN fixed drop number 
case (Table 3, dLWCF and dSWCF) is significantly larger. 
This may be due to the dominance of homogeneous nucleation 
in the base case (and lack of heterogeneous nucleation). 
The results indicate that in both cases the increase in homo- 
geneous nucleation increases the gross cloud forcing, but the 
LW dominates over the SW. There is no evidence of signifi- 
cantly different cirrus AIE with the different ice nucleation 
schemes. Results from ECHAM5-FixedDN simulations are 
similar, but as noted arise from a different response (increasing 
heterogeneous nucleation), with reductions in SW cloud 
forcing resulting in a similar effect. 

[ 40 ] We can also look at how various CAM5 simula- 
tions differ from the CAM5-FixedIN runs. As detailed in 
Section 2.3, differences are denoted by ‘A’ in Table 5, and 
represent differences between the 2000-1850 changes (‘(f) 
in Table 3 between simulations and the FixedIN case. Results 
are presented in Table 5 using 3 different estimates. AIE 
estimates as a difference from the Fixed IN case differ 
between methods by 0. 1-0.2 Wm -2 . This is a significant 
fraction of the mean value, but all estimates have the same 


Table 4. Summary of AIE Estimates From Different Fixed Drop 
Number Runs in Wm -2 as the Change Between Runs With 2000 
and 1850 Aerosol Emissions 11 


Run 

rfCF 

dR -dFSC 

CAM5-LP-FixedDN 

0.31 

0.28 

CAM5-BN-FixedDN 

0.19 

0.30 

ECHAM5-FixedDN 

0.39 

0.45 


“Columns indicate change in cloud forcing (rfCF) and change in TOA 
flux exclusive of clear-sky shortwave change due to direct effect of 
aerosols (dR -rfFSC). See text for details. Simulations from Table 1 . 
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Table 5. AIE as Differences (A) From Fixed IN Runs (in Wm 2 ) a 


Run 

ARFP 

ARFP - AFSC 

ACF 

Mean 

CAM5-LP 

0.22 

0.30 

0.22 

0.25 

CAM5-LP-BC 0.1% 

0.22 

0.27 

0.08 

0.19 

CAM5-BN 

0.27 

0.27 

0.08 

0.21 

CAM5-LP-Low 

0.33 

0.44 

0.17 

0.31 

CAM5-BN-Low 

0.33 

0.46 

0.23 

0.34 

CAM5-LP-W10 

0.56 

0.55 

0.52 

0.53 

CAM5-LP-BC 2% 

0.16 

0.17 

-0.11 

0.08 


“Columns indicate change in TOA flux (R) (or the RFP) between 
two runs (ARFP), the change in RFP exclusive of the clear-sky shortwave 
change due to direct effect of aerosols (ARFP - AFSC), the change in 
cloud forcing (ACF) and the mean of these three estimates. See text for 
details. Simulations from Table 1. The last 2 simulations (W10 and BC 
2%) are sensitivity tests that are not included in the average values 
reported since their basic states are very different. 


sign (positive) and similar magnitude. There is not much 
systematic difference in the ice AIE diagnosed from the dif- 
ferent simulations. The CAM5-LP-Low case, with lower ice 
mass appears to be slightly more sensitive to changes in IN 
(mean AIE of 0.31 Wm - ~) than the CAM5-LP case (mean 
AIE of 0.25 Wm -2 ). The fixed drop number simulations are 
performed with the CAM5-LP-Low microphysical balance, 
and so these estimates are consistent with the fixed drop 
number estimates above. The CAM5-BN case has a slightly 
lower average AIE (0.21 Wm -2 ), probably because of the 
larger fraction of homogeneous nucleation and less hetero- 
geneous nucleation makes it less sensitive while the CAM5- 
BN-Low case has much lower ice number concentrations 
(Table 2), more similar to observations at cold temperatures 
(Figure 2). This simulation with lower ice number con- 
centrations has a relatively large change in ice number con- 
centration due to present-day aerosol emissions (Table 3), 
and consequently a slightly larger ice AIE (Table 5), of 0.23- 
0.46 W m -2 depending on method used to estimate it. 

[ 41 ] Because the ice nucleation parameterizations are sen- 
sitive to vertical velocity, we also perform a sensitivity test to 
understand the impact of different assumptions on ice AIE. 
The CAM5-LP-W10 simulation has vertical velocities which 
are generally higher than those observed in cirrus clouds 
(>0.5 ms -1 ). With significantly higher sub-grid vertical 
velocities (CAM5-LP-W10) Nj is higher (Table 2). Cloud 
forcing is strong (Table 2) and the resulting ice AIE is 
stronger (~0.5 Wm -2 , Table 5). This case is an extreme, and 
illustrates moderate sensitivity to sub-grid vertical velocity, 
but not enough to alter the general conclusions of the paper. 

[ 42 ] Finally we consider the impact of BC on AIE for cirrus 
clouds. To do so, we look at the differences between the 
estimates of AIE derived using variations from the Fixed IN 
runs in Table 5. For the higher efficiency BC cases (CAM5- 
LP-BC 2% and 100%), the mean state is significantly dif- 
ferent from the CAM5-Fixed IN case, so we do not analyze 
the ice AIE quantitatively as a difference from the CAM5- 
Fixed IN run. The mean of three estimates from the CAM5- 
LP case is 0.25 Wm -2 , while that of the CAM5-LP-BC 
0.1% case is 0.19 Wm -2 . The increase in ice number with 
anthropogenic emissions in CAM5-LP of 19 L -1 (Table 3) 
is slightly reduced in CAM5-LP-BC 0.1% (14 L -1 ) due to 
slightly more active heterogeneous BC nuclei. The differ- 
ence in AIE (—0.06 Wm -2 ) is not significant given the 
spread of the different estimates (standard deviation of all 


estimates is 0.10 Wm -2 ). This is the impact of anthropo- 
genic BC in ice nucleation, given an assumed efficiency in 
the LP2007 scheme that maintains a reasonable climate. 
For the CAM5-LP-BC 2% case, the same calculation yields 
a BC effect of —0.17 W m -2 . The impact of BC is of 
the same sign but lower than estimated by Penner et al. 
[2009] with an off-line model, likely due to lower efficiency 
for BC assumed here. 

5. Discussion 

[ 43 ] The simulations provide several different methods 
of looking at the effects of ice nucleation in the current 
climate, and the anthropogenic perturbation to the climate. 
The implementation of the LP2007 and BN2009 scheme in 
CAM5 performed here yield very different results for the 
microphysical balance of clouds. Comparison to in-situ 
observations of ice crystal number and size indicates that most 
of the simulations have slightly higher ice numbers at low 
temperatures. This appears to be due to the influence of 
homogeneous freezing at cold temperatures. We have also 
examined perturbations to the ice nucleation scheme that 
result in very different numbers and sizes, by altering the 
aerosol sizes used in nucleation (CAM5-BN-Low), the overall 
ice mass (CAM5-LP-Low), or the sub-grid vertical velocity 
(CAM5-LP-W10). 

[ 44 ] The parameterizations are constrained with the restric- 
tion that the cloud radiative forcing be similar to observations 
(Section 3). Given this constraint, the BN2009 scheme 
(CAM5-BN) has higher ice numbers and a dominance of 
homogeneous nucleation over heterogeneous nucleation. Two 
different versions of CAM5 with the LP2007 scheme were 
analyzed (CAM5-LP and CAM5-LP-Low), and an alternative 
version with significantly less ice mass (CAM5-LP-Low) 
also has a significant heterogeneous nucleation fraction. The 
balance of homogeneous and heterogeneous nucleation is 
different in ECHAM5-HAM, with much lower ice water 
path. The simulations with lower present-day ice water path 
(ECHAM5 -F ixedDN, CAM5-LP-Low, Table 2) or low ice 
number (CAM5-BN-Low) appear to have slightly higher ice 
AIE (Tables 4 and 5). This indicates that lower mass or 
number implies higher sensitivity to changes. 

[ 45 ] We have also investigated the impact of Black Carbon 
aerosols on ice nucleation in CAM5 (CAM5-LP-BC). Results 
depend entirely on the assumed efficiency of BC as an ice 
nucleus. This is highly uncertain, but most studies show little 
potential for significant heterogeneous nucleation on signifi- 
cant fractions of BC aerosols below the homogeneous 
nucleation limit. This study looked at the impact of assuming 
that different fractions of the BC number were efficient 
heterogeneous ice nuclei. For values of 0.1%, there is a 
small cooling due to reduced increase in N„ consistent with 
other work [Penner et al., 2009], For higher BC efficiency 
(2% or 100%) gross cloud forcing was significantly larger 
than observations (Table 2). Thus a high efficiency is not 
consistent with the CAM5 formulations used here. The BC 
treatment using BN2009 did not impact climate or ice cloud 
microphysics significantly (not shown). 

[46] In all CAM5 cases examined, the ice nucleation balance 
shifts to more homogeneous nucleation as sulfate increases 
from pre-industrial (1850) to present (2000). Seifert et al. 
[2004] found using in-situ data that increases in (sulfate) 
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aerosol number density corresponded to increased ice num- 
ber concentrations. The effect is complex, however this lends 
observational support that an increase in ice crystal number 
from increased sulfate is plausible. The opposite effect was 
noted in ECHAM5-HAM (see below). 

[ 47 ] Using the sensitivity tests as a measure of the uncer- 
tainty, we can summarize the ice AIE results. If all experi- 
ments, models and estimates are treated equally from Tables 4 
and 5 this yields 2 1 estimates from simulations with reason- 
able basic states with a mean of 0.27 ±0.10 Wnf 2 , where the 
uncertainty is one standard deviation. The CAM5-LP-W10, 
CAM5-LP-BC 2% and CAM5-LP-BC 100% experiments are 
not included in this estimate, since their basic states are quite 
different than observations, or other simulations. The estimate 
represents about 20% of the total simulated change in short- 
wave cloud forcing (shortwave AIE) for liquid and ice 
clouds in CAM5 of — 1 .6 Wnf 2 , and is of the opposite sign. 
The indirect effect of cirrus in CAM5 simulations is to drive 
clouds toward increased homogeneous nucleation, 20% higher 
ice number concentration, and an increase in gross cirrus 
cloud radiative forcing, with the long-wave slightly larger 
than the shortwave effect, resulting in a warming. The CAM5 
experiments with more heterogeneous nucleation (CAM5- 
LP-Low) are more sensitive to changing emissions, and have 
slightly higher AIE than the base (CAM5-LP) case, or 
the CAM5-BN case with mostly homogeneous nucleation. 
The effect of adding BC as heterogeneous nuclei slightly 
reduces the warming (—0.06 Wm* 2 ), but not significantly. 

[48] The ECHAM5-HAM results indicate similar AIE 
but for different reasons. In this simulation the ice clouds 
are thinner, with significantly less ice, and the change from 
pre-industrial to present-day is an increase in heterogeneous 
nucleation due to dust coated with sulfate and soluble organics, 
while CAM5 has an increase in homogeneous freezing due to 
sulfate. An increase in heterogeneous freezing can lead to 
more or less ice crystals depending on updraft velocity 
[Karcher and Lohmann, 2003]. In these ECHAM5-HAM 
simulations heterogeneous freezing causes an increase in 
ice crystal number. This nevertheless results in a smaller 
shortwave cloud forcing because the cloud cover decreases. 
This adds to the longwave warming. Thus the increase in 
crystal number makes the LW positive and SW negative, but 
the change in cloud cover has the opposite effect. The com- 
bination for the thin ECHAM5-HAM clouds (low IWP) is 
positive. That the overall results are similar between CAM5- 
FixedDN and ECHAM5-FixedDN is function of the different 
ice water path in the ECHAM5-FixedDN and CAM5- 
FixedDN simulations and the ECHAM5-FixedDN cloud 
cover change. 

[ 49 ] These simulations do not provide complete repre- 
sentation of all ice nucleation mechanisms. Soluble organics 
have been found to be IN when they become solid at 
very cold temperatures [Murray et al., 2010], which might 
help explain low ice crystal numbers and impede ice crystal 
formation [Zobrist et al, 2008]. There are also suggestions 
that oxalic acid [Zobrist et al. , 2006] or crystalline ammonium 
sulfate [Abbatt et al., 2006] could act as heterogenous IN, 
reducing ice crystal numbers. Hendricks et al. [2011] and 
Wang et al. [2010] also found that including BC and other 
heterogenous ice nuclei reduced TV,-. Note that our increases in 
ice crystal numbers are occurring with the increase in sulfate 
from pre-industrial to present-day. Simulations with more 


heterogenous nucleation (e.g, CAM5-LP or preindustrial 
cases) tend to have lower ice crystal numbers, consistent with 
Hendricks et al. [2011], We have examined simulations 
with a wide base state of ice clouds, with varying balances 
of heterogeneous and homogenous nucleation (and ice 
number concentrations), and found consistent results even 
with different mechanisms, lending some confidence that 
‘missing’ present-day IN that altered (lowered) simulated TV,- 
would not significantly alter the results that depend on 
changing sulfate and a more consistent shift toward increased 
crystal number due to homogeneous freezing (CAM5) or 
heterogeneous freezing (ECHAM5-HAM). 

6. Conclusions 

[ 50 ] This work uses 2 GCMs with ice supersaturation and 
ice nucleation, and 3 commonly used ice nucleation schemes 
in large scale models. The GCM simulations generally are 
on the high end of observations of ice crystal number 
concentration at low (tropical cirrus) temperatures. Several 
sensitivity tests have lower ice concentrations consistent 
with observations, but the ice AIE results are similar. 
Different ice microphysical states occur with different 
parameterizations of the ice cloud nucleation process. 
At reasonable BC nucleation efficiencies (0.1%) that are 
consistent with laboratory measurements, BC reduces A) and 
cools slightly. High BC efficiencies are not consistent with 
observed global radiative balance of clouds in CAM5-LP. 
In CAM, indirect effects of cirrus clouds occur mostly due 
to increases in ice crystal concentrations due to homoge- 
neous nucleation from anthropogenic sulfur emissions. 
The resulting ice indirect effects are slightly larger for those 
states with less homogeneous nucleation and smaller ice 
numbers. There is a moderate sensitivity to sub-grid vertical 
velocities assumed for ice (higher ice AIE with higher sub- 
grid velocity and larger TV,). In simulations with low ice 
numbers (CAM5-BN-Low), ice AIE is 20-50% larger, 
suggesting that low ice number states may be more sensitive 
to perturbations. 

[ 51 ] The ECHAM model has similar ice AIE, but for 
different reasons. In ECHAM5-HAM, thinner ice clouds 
and the treatment of immersion freezing of mineral dust 
result in an increase in heterogeneous nucleation and higher 
crystal numbers in present-day because heterogeneous 
freezing occurs at lower ice supersaturations and increased 
sulfate acts to coat dust. In ECHAM5-HAM with thin ice 
clouds (but high crystal numbers) the increase in crystal 
number is outweighed by a decrease in cloud coverage 
causing a reduction in SW cloud forcing (net warming). 
CAM5 can reproduce the positive SW cloud forcing effect 
for the 100% BC case (Table 3) where cloud fraction also 
decreases by over 1%. The mechanisms between the 
two models are different as a consequence of different IN 
parameterizations, but their effect is the same due to different 
mean states (ice water paths) of cirrus clouds. 

[ 52 ] The total ice AIE is estimated at 0.27 ± 0.10 Wm* 2 
(lcr uncertainty ). This represents about 20% of the simulated 
shortwave Aerosol Indirect Effects of —1.6 Wm* 2 for 
CAM5-LP. We note that the simulations neglect possible 
effects of aerosols on convective clouds [e.g., Lohmann, 
2008] that might alter detrained ice crystals. 
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[53] This work highlights that ice indirect effects are 
dependent on the assumed mechanisms for ice nucleation and 
the base state of ice clouds. Regardless of ice microphysical 
balance, addition of anthropogenic emissions, largely of 
sulfate, acts to increase ice number and results in positive 
forcing of consistent magnitude. To better constrain AIE on 
cirrus clouds, better understanding of the role of soot and 
sulfate in heterogeneous and homogeneous nucleation is 
necessary. In addition, better observations of the mean state 
of cirrus cloud microphysics, including ice water path and 
crystal size and number, would help discriminate between 
these realizations in CAM5 and ECHAM5-HAM and better 
constrain the AIE of cirrus clouds. 
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